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FOREWORD 


This publication is a documentation of a wind tunnel test that took place 
in the Naval Surface Weapons Center's Hypervelocity Wind Tunnel #9 in December 
1979. ‘The experimental program was a heat transfer test made at Mach 14 on a 
sphere-cone body instrumented with co-axial surface thermocouples. This test 
was the "trial run" for the use of these gages in the hypervelocity wind tunnel. 


This publication describes the thermocouples used, together with a description 
of how heat transfer rates are calculated from the surface temperature measurements. 
It explains the details of the test set-up, the model configuration, and the data 
reduction technique. It also gives the final results of this test and states 
the accuracy and advantages of this method. 


Special acknowledgements are extended to the Arnold Engineering Development 
Center for their assistance in sending reports that described their experiences 
with the use of co-axial thermocouples in their wind tunnels. The reports helped 
us to avoid unnecesary problems with the implementation of the technique. 
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INTRODUCTION 


. ‘In the Naval Surface Weapons Center's Hypervelocity Wind Tunnel #9, heat 
transfer measurements were generally made using Gardon gages (Reference 1). 
However, the use of co-axial thermocouples to measure heat transfer offers some 
important advantages: | 


1. Calibration stability 

2. Sturdy design 

3. Quick response time . 

4. Ability to be contoured to model surface 


Since these surface thermocouples had never been used in Tunnel 9, a 
shakedown test plan was established to “iron out” any problems associated with 
the use of these gages. This publication is a documentation of that shakedown 
test and its results... _ 

\ 
WIND TUNNEL FACILITY 


The shakedown test was conducted in the Hypervelocity Wind Tunnel #9 from 
10-12 December 1979. Tunnel 9 has a five foot diameter test cell that uses 
nitrogen as the working fluid. For the shakedown test the Mach 14 nozzle was used 
to expand the nitrogen. The average run time for Tunnel 9 is 1.3 seconds, with 
uniform flow occurring during the last 0.7 seconds of the run (Fig. 1). During 
this uniform flow, the model can be pitched through a range of angles of attack. 
More information about Tunnel 9 can be found in Reference 2. 


T ? 
Gardon, Robert, "An Instrument for the Direct Measurement of Intense Thermal 
Radiation," The Review of Scientific Instruments, Vol. 24, No. 5, May 1953. 


2Hill, J. A. F., Wardlaw, A. B., Jr., Pronchick, S. W., and Holmes, J. E., 


"Verification Tests in the Mach 14 Nozzle of the Hypervelocity Tunnel at NSWC 
(White Oak)," AIAA Paper 77-150, Jan 1977. 


QV 
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DESCRIPTION OF CO-AXIAL SURFACE THERMOCOUPLE 


The model number TCS-101-E thermocouples used in the shakedown test are 
manufactured by Medtherm Corporation in Huntsville, Alabama. Figure 2(a) shows 
a picture of a typical thermocouple, and Figure 2(b) shows a cross-sectional 
view of the sensing probe of the thermocouple. The sensing probe consists of 
two metals, chromel and constantain; chromel being the outer tube (first thermo- 
couple element) and constantan being the center wire (second thermocouple element). 
The two elements are insulated except for a vacuum deposited metallic coating 
which is placed on the end of the probe to form a thermal junction between the 
chromel and the constantan. Therefore, temperature readings are measured only at 
the very tip of the sensing probe. Just below the sensing probe is a mounting 
thread so that the thermal junction can be positioned relative to the surface of the 
wall (Fig. 2(c)). 


Surfaces of models tested in the wind tunnel are often curved. Since 
any surface discrepancies could cause disturbances in the boundary layer, the 
thermocouples are contoured to the surface by sanding down the tip of the thermo- 
couple using 180 grit sandpaper. Although this sanding process takes away the 
plating, a thermal junction is still created by the blending of the two metals. 


The average time response for a thermocouple with the vacuum deposited coating 
is one microsecond. For the "blended metal" thermal junction the average time 
response is about ten microseconds. 


Thermal properties for chromel and constantan are given in Table 1, and 
thermal properties for the chromel-constantan thermocouple are given in Table 2. 


THEORY BEHIND CALCULATING HEAT FLUX FROM SURFACE THERMOCOUPLES 


For a one-dimensional heat flux ir.to a homogeneo:'s, semi~infinite solid, 
the heat flux, q(t) can be calculated from the change in surface temperature, T(t), 
from t = 0 by the following equation (References 3 and 4): 


t 
a(t) = K(nk)7? se + T(t) - T(t) | (1) 
t 


where T is the dummy time variable of integration. Since a linear relationship 
is assumed to exist between the actual thermocouple output voltage, E(t), and 
temperature, (AE = 6AT), Equation (1) can be rewritten as: 


Carslaw, H. S. and Jaeger, J. C., Conduction of Heat in Solids, Second Edition, 
Oxford, Clarendon Press, 1959. 


4. ; : 
Vidal, R. J., “Model Instrumentation Techniques for Heat Transfer and Force 


Measurements in a Hypersonic Shock Tunnel," CAL Report No. AD-917-A-1, Feb 1956, 
WADC TN 56-315, AD 97238. 
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ay 1] E(t) 1 0 7° E(t) - E(t) 
q(t) = K(mk) 7 6 oe Hf a a (2) 
t 2 re) (eaqyl* 


Since the integral in Equation (2) is very difficult to evaluate, a method 
will be illustrated later in this report (See DATA REDUCTION) for the calculation 
of q(t). 


PARAMETERS FOR CREATING A HOMOGENEOUS, SEMI-INFINITE SOLID 


Since Equations (1) and (2) are based on the fact that the heat is flowing 
into a homogeneous, semi-infinite solid, there are three parameters to consider 
in making a wind tunnel model wall with a thermocouple mounted into it behave 
as a homogeneous, semi~infinite wall. 


The first parameter to consider is the lumped thermal property, Yk/i, of 
the chromel, constantan, and the model wall. If this property is relatively the 
same for all three materials, then the concept of homogeneity, is valid. Since 
k/K for chromel and censtantan is approximately 2.45 ft*-sec’-°F/BTU, then the 
thermocouple itself is essentially homogeneous. To prevent any radial heat 
conduction the material that the thermocouple is mounted in (model wall) should 
also have aqjJk/K value approximately equal to 2.45. 


The second and third parameters are the duration of the actual wind tunnel 
run and the effective length of the thermocouple sensing probe. If the wind tunnel 
run is of short duration and the sensing probe is long enough, then the semi- 
infinite assumption is valid (Reference 5). Since the duration of an average 
Tunnel 9 run is 1.3 seconds, an appropriate sensing probe length, L, can be 
selected using the graph shown in Figure 3. Therefore, 


L(kt)7? = 2.6 for 0% error. 


8.84 X 10? ta’ /gee (constantan) 
t 1.3 seconds 


then L > .28 inches. 


The appropriate sensing probe length and wall thickness for a model in Tunnel 9 
should be greater than .28 inches. 


MODEL CONFIGURATION 
The model configuration that was tested was a sphere-cone type body. The 
nosetip had a 1.8" radius and the cone half-angle was 7. The model was made 


out of 17-4°H Stainless steel. The Yk/K for this material is 2.44 ft2sec -°F/BTU,* 
which is very close to theqk/K of the chronel-constantan thermocouple. 


Brown, H. K., "The Theoretical Response of Heat Transfer Gages Employed in Shock 
Tubes,"' AVCO Research Laboratory, Research Note 58, Feb 1958. 


* e 
Obtained from Materials Selector 75, Vol. 80, No. 4. 
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There were two interchangeable conical sections used in the test. The first 
conical section was referred to as the “thick wall body" because its wall was 
3/8" thick (which is thicker than the critical 0.28 inches), and the second conical 
section was referred to as the "thin wall body" because its wall thickness was 
only 0.125" thick. Figure 4 shows a sketch of the two configurations. 


INSTRUMENTATION 


In the nosetip of the model, two co-axial thermocouples were mounted as shown 
in Figure 5. Thermocouple "1" was mounted directly in the wall; the wall at 
that point being thicker than 0.28 inches. However, co-axial thermocouple "2" was 
mounted in the wall inside a 17-4PH stainless steel 0.5" diameter plug that was 
required to make the wall thicker than 0.28 inches. 


In the "thick wall" conical section, three co-axial thermocouples and three 
Gardon gages were mounted as shown in Figure 6. The three thermocouples were 
mounted 5.83 inches downstream from the nosetip; one thermocouple on the leeward 
meridian, one on the 90° meridian, and one on the windward meridian. Each Gardon 
gage was mounted one inch downstream from the thermocouples; one on each of the 
meridians. 


In the "thin wall" conical section, three co-axial thermocouples and three 
Gardon gages were also mounted in the same positions as the "thick wall" body, 
as shown in Figure 6. However, since the wall was only 0.125" thick, the thermo- 
couples were mounted in the wall with plugs that would make the wall 0.375". The 
plugs had varying diameters to determine a minimum permissible plug diameter, 


The Gardon gages used in both the "thick" and “thin" wall bodies were 
manufactured by Thermogage and had been used in previous wind tunnel tests. 
Each gage's heat flux sensitivity, C, was calculated using a calibrated lamp as 
a known heat source. Each gage's time delay constant, Tg, used in the data 
reduction equations (see DATA REDUCTION) was then calculated by observing the time 
it took for each gage to respond to 63.2% of its fullscale output for a step heat 
input. The Gardon gages were used in this shakedown test as a check to the co-~ 
axial thermocouples. 


TEST SCHEDULE 


The test matrix and run conditions are given in Table 3. The pitch sweeps 
were set up to compare upsweep (Run 496) with downsweep (Run 498) data, to compare 
a static angle of attack (Run 497) with the upsweep and downsweep data, and to 
compare thick and thin wall configurations (Runs 496 and 499). 


DATA REDUCTION 


As was stated previously, Equation (2) is very difficult to evaluate. For 
reduction of the raw surface thermocouple output, E(t), into heat flux data, the 
Dixon Method (Reference 6) was used. The Dixon Method is a two-step procedure 


6 ; ‘ 

Kendall, D. N. and Dixon W. P., "Heat Transfer Measurements in a Hot Shot Wind 
Tunnel," presented at the IEEE Aerospace Systems Conference, Seattle, Washington, 
11-15 Jul 1966. 
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that does not require any initial smoothing of the raw thermocouple output. First, 
the total heat transfer to the surface is calculated using the following equation: 


n E(t._,) + E(t.) 
Q(t) = K(km) 267! ——-_—_, At (4) 
; ( 


i=l it = t._)) + (to = t.) 
where n= 0, 1, 2... (t/At + 1) and where At is an equal time increment. 


Then, the heat transfer rate is computed by differentiating Q(t): 


q(t) = ar 


The expression for differentiating Q(t) is described in Reference 7 and is: 


dQ(t) 1 
a(t) = —ge— = TOyTHEY [2a ~ QE Ag) * UC) * 2) v 


A sample voltage was recorded just prior to the wind tunnel run. This 
sample voltage was then subtracted from all subsequent voltage readings. Therefore, 
at t, = 0, E(t) = 0 which imples that q(t.) = 0. 


For the reduction of the Gardon gage output, E(t), the raw data was first 
smoothed, reversed, and smoothed again using a sixth order Butterworth digital 
filter set at a cutoff frequency of 5Hz. Heat transfer rates were than calculated 
using the following standard Tunnel 9 equation: 


: = dE(t) 
q(t) =C [ece) + 1% He) ] (7) 
where c = calibrated gage sensitivity Ctey? 
T> = calibrated time delay constant 
Th dE(t) . . . 
e term —-— is calculated by the method given in Reference 7. 
dE(t) 1 
8 aie [-28¢t ~ E(t _,) + E(t 4,) + 28(Ea49) (8) 


where n= 0, 1, 2 2 Cat + 1) 


‘enrich, Fredric F., "Differentation of Experimental Data Using Least Squares 
Fitting," Journal of the Aeronautical Sciences, Vol. 22, No. 2, Feb 1955. 
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Equation (7) is only valid if at t = 0, q(t.) = 0. Therefore, a sample of 
data was recorded just prior to each wind tunnel run, and this sample voltage 
was then subtracted from all subsequent voltages.so that at t_ = 0, E(t _) = 0 
implying that q(t) = 0. . : 


From the heat transfer rates calculated from the co-axial thermocouple and 
Gardon gage readings, Stanton numbers were calculated by the following equation: 


ST = 4 [Potato - 7 | os (9) 


where q = calculated heat transfer rate (BTU/£t?-sec) 


0,, = free stream density (1bm/ ft?) 


[=] 
" 


free stream velocity (ft/sec) 


BIU fe) 


/lbm - FP 


Q 
i] 


heat capacity for nitrogen = 0.2481 


F : ey 
Ty, = equivalent ideal gas supply temperature (“F) 
(calculated from T, and tables in Reference 8) 


* 
= measured wall temperature (°F) 


=) 
i 


The free stream properties are calculated from a pitot tube measurement in 
the flow and a supply pressure, Ps? measurement. 


RESULTS 


Table 4 is a listing of the data obtained from the shakedown test. It should 
be noted that T5 went bad on Run 499, and G3 went bad on Run 498. The listing 
only shows data during the "uniform flow" portion of each run. Figures 7 through 
14 show plotted data of Stanton number vs. angle of attack for all four runs. 

Heat transfer calculations made by the G.E. 3-D Viscous Code (Reference 9) are 
also shown on these figures. 


B : : : 
Cullotta, S. and Richards, B. E., "Methods for Determining Conditions in Real 
Nitrogen Expanding Flows," VKI-TN-58, Feb 1970. 


* 
For Gardon gage data the nearest co-axial thermocouple temperature reading was 
used as the t value. 


Hecht, A. M., Nestler, D. E., and Richbourg, D. H., “Application of a Three- 


Dimensional Viscous Computer Code to Reentry Vehicle Design,"' AIAA Paper 79-0306, 
Jan 1979. 
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ACCURACY 


F Comparisons will be made with respect to the repeatability of the upsweep 
(Run 496) and downsweep (Run 498) data; the repeatability of the upsweep (Run 496), 
downsweep (Run 498), and static angle (Run 497) data; and the repeatability of 

the "thick wall" configuration (Run 496) and the "thin wall" configuration (Run 
499) data. A comparison will also be made between the calculations of the G.E. 

3-D Viscous Code (Reference 9) and the data for Runs 496, 498, and 499. It should 
be noted that the following tunnel properties have the following previously 
observed accuracies: 


Supply pressure, PS - +.4% 


Supply temperature, T> - -1.7% to +.5% 

Pitot measurement - +.3% 

Free stream Mach number, M, - +.4% 

Free stream pressure, P, - +2.8% 

Free stream unit Reynolds number, Re, /ft - -1.4% to +2.8% 


The angle of attack measurements are accurate to within 0.1° for Run 496 
and to within 0.3° for Runs 497 through 499. 


The Dixon method calculates heat transfer rates to within an accuracy of 
less than 1%. As for the Gardon gages, the gage sensitivities and the time delay 
constants are accurate to +54. 


Comparison of Upsweep vs. Downsweep Data. Table 5 lists the accuracies for 
the repeatability of the Stanton number data for Run 496 (upsweep) vs. Run 498 
(downsweep) for each gage at 5 angles of attack. The repeatability for the two 
runs shows an average percentage difference of about 7.84. 


Comparison of Upsweep~Downsweep vs. Static Angle Data. Table 6 lists the 
accuracies for the repeatability of the Stanton number data for Run 496 (upsweep) 
vs. Run 497 (static angle) and Run 498 (downsweep) vs. Run 497 (static angle) 
for each gage at an angle of attack of 10°. The repeatability between the dynamic 
and static data has an average difference of about 4.1%. 


Comparison of "Thick Wall" vs. "Thin Wall" Data. Table 7 lists the accuracies 
for the repeatability of the Stanton number data for Run 496 (thick wall) vs. Run 
499 (thin wall) for each gage at five angles of attack. G2 is not listed because 
it was slightly recessed in the model wall and was, therefore, measuring lower 
heating rates. 13, T4, and T5 were mounted in plugs of 3/4", 1/2", and 1/4" 
diameters, respectively. Each plug mad: the wall thickness 3/8". Since T3 went 
bad, the results are inconclusive as to the minimum diameter plug that can be 
used so that the wall will be semi-infinite in the radial direction. A plug may 


" See footnote 9 on page 10. 
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not be needed as long as the sensor length is greater than 0.28 inch. However, 

T2 was mounted in a 1/2" diameter plug and its repeatability difference throughout 
the test was about 6.5%. Therefore, a configuration of a 1/2" diameter plug and 
thermocouple is a possible working configuration. 


Comparison of Runs 496, 498, and 499 vs. the G.E. 3-D Viscous Code. Table 8 


lists the accuracies for agreement of the Stanton number data between Runs 496, 

498, and 499, and the GE 3-D Viscous Code (Reference 9) for each gage at 0 and 

5 angle of attack. The average difference in agreement with the code is about 

8.8%, with the code's calculation of the leeward heating contributing to most of 
the error. 


CONCLUSIONS 


In comparing the surface thermocouple method of measuring heat transfer to the 
use of Gardon gages in Tunnel 9, the thermocouples have distinct advantages: 


1. Gardon gages require a calibration, whereas the thermocouples have an 
inherent bi-metallic calibration. 


2. Gardon gages have a slow response time (on che order of 50 msec) that 
must be rectified in the data reduction procedure to acquire accurate timewise 
data, whereas the thermocouples have an almost "instantaneous" response time 
(about 10 usec). 


3. Gardon gages have a delicate, thin skin that can be broken by the flow 
(e.g., G3 on Run 498), whereas thermocouples are a solid piece of metal that 
cannot be disturbed by the flow. 


4. Gardon gages cannot be contoured to the model surface, e.g., G2 was 
slightly recessed in the model wall causing it to measure a lower heating rate, 
whereas the thermocouples can be contoured exactly to the model surface. 


An estimate of the accuracy of heat transfer rates by the thermocouples 
is +6%. This is slightly better than the 7% accuracy that has been observed for 
Gardon gages in Tunnel 9. In light of this accuracy along with the advantages 
over the Gardon gages, the co-axial thermocouples proved to be a viable method 


for measuring aerodynamic heating during pitch sweeps in Tunnel 9. 


see footnote 9 on page 10. 
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FIGURE 2(a) A TYPICAL TCS MODEL THERMOCOUPLE 


FIRST THERMOCOUPLE 
ELEMENT (TUBE: 


INSULATION 


SECOND THERMOCOUPLE 
ELEMENT (a WIRE, 


vVaCUUM OEPOSITED PLATING 


FIGURE 2(b) CROSS - SECTIONAL VIEW OF FIGURE 2(c) THERMOCOUPLE 
THERMOCOUPLE. (PLATING 


PROBE MOUNTED 
AND INSULATION EXAGGERATED IN MODEL WALL 
IN SIZE) 
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THIN WALL BODY 


FIGURE 4 MODEL CONFIGURATIONS 


16 


NOILVINAWNNYLSNIdILaSON S 3HNDIS 
31V9S O1 LON 


M3IA LNOYS 


M3IA 3QGis 


(GYVMQGNIM) .0 


NSWC MP 80-151 
17 


Y¥L8°L 


NSWC MP 80-151 


THICK WALL SECTION 180° 


FROM NOSETIP 
5.83” 0° (WINDWARD) 


THIN WALL SECTION 180° 


0° (WINDWARD) 


FIGURE 6 CONICAL SECTION INSTRUMENTATION (NOT TO SCALE) 
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* kk 
TABLE 1 THERMAL PROPERTIES OF CHROMEL AND CONSTANTAN 


Thermal Conductivity 3 


K (BTU/in-sec-°F) @ 75°F .242 X 10 3 


.267 X 10 


Specific Heat 
C. (BTU/1b-°F) @ 68°F 107 094 


Density 


6 (1b/in?) 


315 .322 


Thermal, Diffusivity 


k (in*/sec) 7.18 X 10> 8.84 X 10> 


Melting Point 
(°F) 2600 2210 


* 
Properties given by the Hoskins Manufacturing Company, Detroit, Michigan 48208 


kk 
Properties given by the Thermo-Electric Company, Inc., Saddle Brook, 
New Jersey 07662 


NSWC MP 80-151 


* 
TABLE 2 THERMAL PROPERTIES OF CHROMEL-CONSTANTAN THERMOCOUPLE 


Thermal Property Value 
Lumped Thermal Property 
{kK/K (£t2-sec + °F/BTU) 
. . . : ** 
Thermoelectric Sensitivity 
& (uv/°F) 


Heat Flux Sensitivity 


5 
none (alse —) 
BIU/ft© - sec 


¥ 
Theoretical values 


*k 
Value obtained from National Bureau of Standards Circular #561 for a chromel- 
constantan thermocouple 


TABLE 3 TEST SCHEDULE AND RUN CONDITIONS 


Angle of Attack 
Model Average Average Average Average Pitch Sweep dur- 
Configuration Mach No. i To (°F) Re,/ ft ing Uniform Flow 
"Thick Wall" 19604 3.44% 10° 0° for .4sec then 
Body a pitch from 0 
to 18 
"Thick Wall" 19538 3.73 X 10° =10° constant 
Body 


6 


Pitch from 19° to 
Body 22° to 0 


"Thick Wall" 19782 3.92 X 10 


"Thin Wall" 20741 3.85 X 10° Same as 496 


Body 
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TABLE 5 ACCURACIES FOR REPEATABILITY OF RUN 496 VS. RUN 498 
(UPSWEEP VS. DOWNSWEEP) 


(Note: Values in % difference in agreement) 


TABLE 6 ACCURACIES FOR REPEATABILITY OF RUN 496 AND RUN 498 VS. RUN 497 
(DYNAMIC SWEEP VS. STATIC) 


3.8% 19.3% 5.5% 15.1% 


(Note: Values in % difference in agreement) 
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s 


TABLE 7 ACCURACIES FOR REPEATABILITY OF RUN 496 AND RUN 499 
("THICK WALL" VS. "THIN WALL") 


* 
G2 was recessed in model wall on Run 499. No comparison of this gage was made. 


(Note: Values in % difference in agreement) 


TABLE 8 ACCURACIES FOR AGREEMENT OF RUNS 496, 497, AND 499 vs. 
THE G.E. 3-D VISCOUS CODE 


(Note: Values in % difference in agreement) 
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TERMS 
angle of attack (°) 
calibrated Gardon gage sensitivity (4) 
specific heat 
output voltage 
thermal diffusivity 
thermal conductivity 
minimum effective sensing probe length 
free stream Mach number 
free stream Mach number 
supply pressure (psia) 
free stream pressure (psia) 
supply pressure (psia) 
heat transfer rate (BTU/£t?-sec) 
cumulative heat transfer to a surface (BTU/£t?) 
heat transfer rate (BTU/£t?-sec) 
free stream Reynolds number 
free stream Reynolds number 
free stream density (1bm/£t>) 
Stanton number 


time 
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NSWC MP 80-151 


TERMS (Cont. ) 
temperature 
supply temperature (°F) 
equivalent ideal gas supply temperature (°F) 
measured wall temperature (°F) 
free stream temperature (°F) 
supply temperature (°F) 
equivalent ideal gas supply temperature (°F) 
measured wall temperature (°F) 
free stream velocity (ft/sec) 
free stream velocity (ft/sec) 


angle of attack 


thermoelectric sensitivity C705) 


density 


Ibm ¢,3) 


free stream density ( 
dummy variable of integration 


calibrated Gardon gage time delay constant 


